QUASI-PERIODIC AND PERIODIC SOLUTIONS OF THE TODA 
LATTICE VIA THE HYPERELLIPTIC SIGMA FUNCTION 



YUJI KODAMA, SHIGEKI MATSUTANI, EMMA PREVIATO 



Abstract. M. Toda in 1967 (J. Phys. Soc. Japan, 22 and 23) considered a lattice 
model with exponential interaction and proved, as suggested by the Fermi-Pasta-Ulam 
experiments in the 1950s, that it has exact periodic and soliton solutions. The Toda 
lattice, as it came to be known, was then extensively studied as one of the completely 
integrable (differential-difference) non-linear equations which admit exact solutions in 
terms of theta functions of hyperelliptic curves. In this paper, we extend Toda's original 
approach to give hyperelliptic solutions of the Toda lattice in terms of hyperelliptic 
Kleinian (sigma) functions for arbitrary genus. The key identities are given by generalized 
addition formulae for the hyperelliptic sigma functions (J.C. Eilbeck et al, J. reine 
angew. Math. 619, 2008). We then show that periodic (in the discrete variable, a 
standard term in the Toda lattice theory) solutions of the Toda lattice correspond to the 
zeros of Kiepert-Brioschi's division polynomials, and note these are related to solutions 
of Poncelet's closure problem. One feature of our solution is that the hyperelliptic curve 
is related in a non-trivial way to the one previously used. 
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1. Introduction 

The Toda lattice is an "algebraically completely integrable system". As such, it ad- 
mits classes of solutions parametrized by Jacobi varieties of compact Riemann surfaces 
(or algebraic curves), the "algebro-geometric solitons" |GHMTj . One advantage of the 
algebro-geometric solution is that the time flows become linear on a hyperelliptic Jaco- 
bian, and the difference operator is translation on the Jacobian. These are the objects of 
concern here. 

On the other hand, recent work [?, IMPlj on Kleinian sigma functions has focused on 
addition formulae on a certain stratification of the Jacobian; given the relevance of this 
stratification in terms of the orbits of the Toda lattice |AMt lAHMt \CK\ IKoj . our program 
is to study the explicit relationship between the two, with the result of an exchange of 
knowledge. 

In this first paper we identify and study the (quasi-)periodic solutions of the Toda lattice 
equation in terms of the hyperelliptic sigma functions. Our approach is somewhat different 
from the one that exists in the literature, and in particular it gives us a different 'spectral 
curve' and algebraic conditions for periodicity. To give a sketch of our method and results, 
we first say a few words about the history and the significance of the cx-function, especially 
in the field of differential/difference equations. In the 19th century, it was the study of 
Riemann- surface theory that led to the associated algebraic (meromorphic) or analytic 
(theta) functions and the (differential) equations that they satisfy, apart some exceptional 
cases In keeping with such philosophy, Baker in |B3] discovered the KdV hierarchy and 
KP equations for every hyperelliptic curve of genus g as an application of Kleinian sigma 
functions (1903) and posed this challenge: 

These equations put a problem: To obtain a theory of differential equations 
which shall shew from them why, if we assume 

PXfi{u) = -d 2 \ogo(u)/du x du^ 

the function o~(u) has the properties which a priori we know it to possess, 
and how far the forms of the equations are essential to these properties. 

Baker extensively studied the differential equations satisfied by the a function. In the late 
1960s and early 1970s, on the contrary, many authors started from the ("completely inte- 
grable" ) differential equations, and arrived at a spectral curve which is a Riemann surface 
and carries a "Baker- Akhiezer" eigenfunction for the (linearizing) "Lax-pair" equations. 
Mumford in his book on theta functions |Mu] demonstrated the close relationship between 

^KP appears in the case when the afhnc model of the hyperelliptic curve has two points at infinity [?] . 
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the differential equations and the algebraic approach, together with the abelian function 
theory of the 19th century based on Jacobi's theory; he gave an explicit dictionary be- 
tween the theta function for hyperelliptic curves and three polynomials in one variable, 
sometimes called the Mumford triplet (whose definition he attributes to Jacobi), which 
parametrize the Jacobi variety, cf. Remark 15.71 

Recently the Kleinian sigma function was reprised, generalized and studied by several 
authors [BEL, EELJ. These authors showed that sigma is more efficient than Riemann's 
theta function (in fact, the sigma function approaches the Schur polynomials in the limit 
when the curve becomes rational) to solve differential equations. This is our point of view. 

We add, as suggested by the referee, a comment on the geometric significance of sigma: 
for more detail we refer to the aforementioned studies. Since sigma vanishes on the 
(higher) Abel images of the curve X g , it is better suited to be expanded in the abelian 
variables which correspond to the hyperosculating flows (Section EJ) to the image of X g , 
which are the flows of the KP hierarchy. It is in these variables that sigma equals a Schur 
polynomial up to higher-order terms [N] . In fact expansions and computer-algebra work 
have enabled guesses and proofs of the addition formulas which generalize (12. 4p and are 
essential to solving integrable equations. Lastly, sigma (unlike Riemann's theta, cf. |Mu[ 
II. 5]) is invariant under the action of the modular group on the period matrix of X g . 

We also comment on the fact that our addition formulae, which are key to the solution as 
explained in the next paragraph, are particular to hyperelliptic curves. Of course a ddition 



theorems hold for general curves, cf. e.g. Theorems 9.1 and 10.1 in |EEMOP2] for the 
trigonal case; however, they will not be as expedient as the hyperelliptic ones, for example 
an expression in terms of algebraic functions on the curve is as yet unwieldy, as perhaps can 
be expected, since the non-hyperelliptic Neumann systems have polynomial Hamiltonians 
of a very complicated type JSTJ ISU [S3]. The basic reason why the hyperelliptic case 
simplifies, and in fact was the one for which Baker was able to study explicitly the sigma 
function, is the hyperelliptic involution which lifts to the Jacobian in such a way that it 
acts on sigma only by a sign. 

In this article, we consider the solutions of the, one- and two- (time-) dimensional, Toda 
one- (space-) dimensional lattice. Toda gave an exact solution for the Toda lattice equation 
using an identity of elliptic-function theory [T]. For an arbitrary hyperelliptic curve of 
genus g, we construct a meromorphic function on the Jacobian of the curve that obeys the 
Toda lattice equation over the Jacobian of the curve, by using an identity of hyperelliptic 
abelian functions (see Theorems 15.41 and 15 . 51 for the one- and two- (time-) dimensional Toda 
lattice equations respectively). Remark 15. 71 shows that the two-(time-)dimensional Toda 
equation is equivalent to a relation which generates Fay's trisecant formula and Baker's 
derivation of the KdV hierarchy and the KP equation for every genus. Despite this 
equivalence, we stress again here the geometric nature of sigma compared to theta. That 
nature gives a dictionary between abelian and algebraic functions of the curve, specifically, 



Schur polynomials were used by Sato to define his tau function and originally construct the universal 
solution to the KP hierarchy. The relationship of sigma with tau is pursued in [EEG . 
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the affine coordinates of its planar representation, as well as algebraic coefficients for the 
relations in the ring of differential operators that act on the Jacobian (more details are 
given in Remark 14.131 and 15.71 (2) (b)). Weierstrass recognised this fact in his study of 
hyperelliptic theta functions |W] : in fact, he defined the al and Al functions (named 
after Abel) which can be used to give solutions of finite-dimensional Hamiltonian systems 

mm- 

We obtain quasi-periodidj (hyperelliptic) solutions of the Toda lattice, since a priori 
the discrete variable has no rational relation with the period lattice of the curve. Finding 
periodic solutions is equivalent to giving a torsion point of the curve, which satisfies a 
"division polynomial", known as Kiepert's or Brioschi's polynomial. Using a zero of the 
division polynomial (if it exists), we construct a hyperelliptic solution of Toda of genus g 
with period N(> g) in Theorems 16.31 and 16.111 We illustrate the g = 1 case as an example 
and in that case point out that the relation between the division polynomial and Toda 
lattice is the same as Poncelet's closure (Appendix). 

Acknowledgements. One of authors (S.M.) thanks Boris Mirman for bringing Poncelet's 
problem to his attention and Akira Ohbuchi who drew his attention to Galois's study on 
fifth elliptic cyclic point. Y.K. is partially supported by NSF grant DMS-0806219. E.P. 
acknowledges very valuable partial research support under grant NSF-DMS-0808708. We 
are indebted to the referee for corrections and suggestions. 

2. Genus-one case 

In this Section, we demonstrate how one gets an elliptic solution of the Toda lattice 
[Tj IKvAJJ IMat3j . Let Xi be an elliptic curve given by 

X\ : y 2 = x 3 + \ 2 x 2 + A]£ + A 

(2.1) = (x - ei)(x - e 2 ){x - e 3 ), 

where the e's are distinct complex numbers and A2 = — (ei + &2 + 63) = 0. 

The Weierstrass elliptic o function associated with the curve X\ is connected with the 
Weierstrass p and ( functions by 

d 2 d 

(2.2) p{u) = - — logo-(u), C{u) = — logo-(u). 

du z du 

The Jacobian of the curve X x is given by Ji = C/ (Zu/ + Zu") using the periods (u/, u"), 
complete integrals of the first kind. The abelian coordinate u, actually defined on the 
universal cover of the Jacobian J\ but used as customary up to congruence when no 
ambiguity arises, is given by 

dx 

^ "=i v 

3 As customary in the literature on the Toda lattice, cf. the monographs [AvMVj and [V]) the word 
"periodic" does not refer to the time, but rather the space variable of the lattice sites. 
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with x(u) = p(u), 2y{u) = p'(u) and oo the infinity point of X\. 

The key to obtain a p-function solution of the Toda lattice is the addition formula, 

a(u + v )a(u — v) 

2.4 p{u) - p(v = - \ > 1 } . 

[a(v )a[u)\ z 

By differentiating the logarithm of (12. 4 j) with respect to u twice, we have 
d 2 

(2.5) - — \og[p(u) - p(v)\ = p(u + v)- 2p(u) + p(u - v). 

Thus for a constant number w and any integer n, by letting u = rax + 1 + to, v = uq, 

(2.6) 
d 2 

~1^2 lo sb( nu o +t + to)- p{u )] = [p((n + 1)« + t + t ) - p{uo)} 



dt 



- 2{p(nuo + t + t ) - p(u Q )] + [p((n - l)it + * + *o) - p(u )]- 



Which, by letting V n {t) := -p{nu + t + t ), V c := -p{u ) and q n := - log[V n (t) - V c ], 
becomes 

d 2 

(2.7) - — q n = e~ qn+1 - 2e~ gn + e" 9 "" 1 (nGZ). 

at 2 

This can be identified with the continuous Toda lattice equation, where q n are the inter- 
action potentials. Indeed, by letting q n = Q n — Q n -i, where Q n is the displacement of the 
n-th particle and obeys the nonlinear differential equation of the exponential lattice [Tj, 

d 2 

(2.8) - — Q n = e Q "- Qn+1 - e ' 1 - 1 " ' 1 (nGZ). 

To provide the connection between Toda's original solution for e\, e2, &3 G R [T] and this 
elliptic solution derived from the addition formula ( 12 Ah with ( 12. 21) . we observe that, by 
letting t = ~ w", 

p(t + nuo - lu") - p(t + nuo) = (ex - e 3 ) ns 2 ((ei - e 3 ) 1/2 (t + nu )) + e 3 , 

where ns(w) = ksn(u + zi^') and the modulus of the Jacobi elliptic functions is k 2 = 
(e2-e 3 )/(ex-e 3 ) [WWl 22-351]. 



3. HYPERELLIPTIC CURVE X g AND SIGMA FUNCTIONS 

In this Section, we give background information on the hyperelliptic 9 functions and 
the a function, a generalization to higher genus of the Weierstrass elliptic a function. 
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3.1. Geometrical setting for hyperelliptic curves. Let X g be a hyperelliptic curve 
defined by 

X g : y 2 = f(x) := x 29+1 + \ 2g x 29 + • • • + A 

where X/s are complex numbers, together with a smooth point oo at infinity. Let the 
affine ring related to X g be R g := C[x, y]/(y 2 — f(x)). We fix a basis of holomorphic 
one-forms 

t x % ~ x dx 



v 



2y 



i = 1, 



We also fix a homology basis for the curve X g so that 

9 9 

H x (X g , Z) = Zaj © Z# 

3=1 0=1 



31 



with the intersections given by [c^, aj] = 0, [f3i, (3j] = and [a^^j] = — aj] = 5ij. We 



consider the half-period matrix u 



U= 2 



UJ 

uj" 



V 



of X g with respect to the given basis where 



to 



Let A be the lattice in C 9 spanned by the column vectors of 2u' and 2u" . The Jacobian 
variety of X g is denoted by J g and is identified with C 9 /A. For a non- negative integer k, 
we define the Abel map from the fc-th symmetric product S k X g of the curve X g to C 9 by 
first choosing any (suitable) path of integrationB 

w : S k X g C 9 , w((x h yi),..., (x k , y k )) 



i=i 



'Vx 



By letting the map k be the natural projection, 



K : C 9 



Jc 



ti- 



the image of k o w is denoted by W4 = n o w(S k X g ). The mapping k o w is surjective 
when > (7 by Abel's theorem, and is injective when k = g if we restrict the map to the 
pre-image of the complement of a specific connected Zariski-closed subset of dimension at 
most g — 2 in Jj g , by Jacobi's theorem. 



The results presented below are independent of the particular choice. 
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3.2. Sigma function and its derivatives. We define differentials of the second kind, 

1 2g-j 

= Y S ( fc + 1 ~ j)^k+i+jX k dx, (j = 1, . . . , g) 



k=j 



and (half of) complete hyperelliptic integrals of the second kind 



V 



v 



For this basis of a 2g-dimensional space of meromorphic differentials, the half-periods 
uj', uj", T]', rj" satisfy the generalized Legendre relation 

-In \ ™VT ^ ( 



(3.1) 



o 



m 1 



o 



uj' uj" 







rj rj 
tics r 

'.z:Ti 



where 9JI = ( ~, ~ „ \. Let T = uj' 1 u" . The theta function on C 9 with "modulus" T 
7 rj J 

and characteristics Ta + b for o, 6 6 C 9 is given by 

1 t, 



X] ex p 



2 tt/ { - t (n + a)T(n + a) + \n + a)(z + 6) 



The a-function ( [Blj . p. 336, |BEL] ). an analytic function on the space C 5 and a theta 
series having modular invariance of a given weight with respect to Wl, is given by the 
formula 



a(u) = 7o exp 



^rj'cu'^u } 9 



uj' 1 u ; T 



where 5' and 6" are half-integer characteristics giving the vector of Riemann constants 
with basepoint at oo and 70 is a certain non-zero constant. The cx-function vanishes 
exactly on /t _1 (W 9 _i) (see for example [Bit XI. 206]). The Kleinian p and ( functions are 
defined by 

Let {(pi(x,y)} be an ordered set of C U {oo}-valued functions over X g defined by 

{x % for i < g, 

x l(i-g)/z\+9 f or i > g anc [ 1 _ g e ven, 
x K i -g)/ 2 \y f or i > g anc i 2 — g odd. 

We note that \4>i(x, y)} give a basis of the (infinite-dimensional) C vector space R g . 

Following |0[ IMP2] . we introduce a multi-index t| n . For n with 1 < n < g, we let \\ n be 
the set of positive integers i such that n + 1 < i < g with z = n + 1 mod 2. Namely, 

n + 1, n + 3, . . . ,g — 1 for (7 — n = mod 2 
n + 1,71 + 3, ...,<7 for g — n = 1 mod 2 
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and partial derivatives over the multi-index tj 




For n > g, we define t| n as empty and cr^ as a itself. The first few examples are given in 
Table 1, where we let jj denote t] 1 and b denote \] 2 . 

Table 1 



genus 


cr,j = 


0"|> = 0"^2 


(71,3 


(7^4 


0^ 5 


0il 6 


oy 




1 


a 


a 


O" 


cr 


(J 


cr 


a 


a 


2 


02 


a 


a 


cr 


O" 


cr 


a 


a 


3 


02 


03 


a 


a 


O" 


(7 


a 


a 


4 


024 


0-3 


(74 


a 


a 


cr 


a 


a 


5 


024 


0"35 


0"4 


05 


a 


a 


a 


a 


6 


0246 


035 


046 


0-5 


06 


a 


a 


a 


7 


0246 


0357 


046 


057 


06 


o- 


a 


a 


8 


02468 


0357 


0468 


057 


068 


a 7 




a 



For u G C 9 , we denote by w' and u" the unique vectors in R 9 such that 

it = 2 W + 2 W. 

We define 

L(u,u) =*u(2V« , + 2VV), 

= exp {2v™(Y 5" - + \H'l") } (g {1, -1}) 

for u, t> G C 9 and for £ (= 2 Vf + 2V'f) G A. Then <r r (u) for u G « _1 (>Vi) satisfies 
the translational relation (|Qj. Lemma 7.3): 

(3.3) a r (w + £) = xM^»(tt)expL(tt+ for w G k -1 (Wi). 

Further for n ^ g, we note that 0^(— w) = (— l) n5,+ 5"( n ~ 1 V[ 1 n('u) for w G K~ l (W n ), 
especially, 

{(Tb(— = — adu) for m G k~ x (W2) 

^(-u) = (-1)% b (m) for u G kT^Wi) 
by Proposition 7.5 in jo] . 

4. The addition formulae 

In this Section, we give the addition formulae of the hyperelliptic a functions which 
are the generalization of the addition formula (12 .4p . These are essential to construct the 
hyperelliptic solution of the Toda lattice. 
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4.1. Generalized Probenius-Stickelberger formula. We first recall the generalized 
Frobenius-Stickelberger formula which gives a generalization of the addition formula (12. 4p . 

Definition 4.1. For a positive integer n > 1 and (xi,yi), . . . , (x n ,y n ) in X g , we define 
the Frobenius-Stickelberger determinant [MPlj . 



0i(^i,2/i) 
0i(^2,2/2) 



(f> n -i{xi,yi) 
4>n-i(x 2 ,y 2 ) 



<f>n-2{Xn-l,Vn-l) 4>n-l{x n -l, Vn-l) 
4>n~2{Xn,y n ) 4>n-l{Xn,y n ) 



1 <f>i(x n -i,y n -i) ■ 

1 4>l{x n ,Vn) 

where <fii(xj,yj) 's are the monomials defined in VS.ty . 
Then we have the following theorem (Theorem 8.2 in [o]): 

Proposition 4.2. For a positive integer n > 1, let (x\, y%), . . . , (x n , y n ) in X g , and 
in k _1 (Wi) be points such that ixW = w((xi,yi)). Then the following re- 



u 



(i) 



u 



(n) 



lotion holds: 



(4.1) 



^(Er=i wW )rii<i^(w 



e n *n((^i,2/i) 



w/iere e n = (-l)9+"(«+i)/2 for n < g and e n = (-l)( 2 "-ff)0?-i)/ 2 forn>g + l. 

4.2. The algebraic addition formula. We first describe linear equivalence of divisors 
on the curve X g (which will result on the addition law on the Jacobian) by algebraic 
formulas [MPT] : 

Definition 4.3. For given Pi, . . . , P n G we define 

^ n +l(P,Ql,...,Qn) 



^ n (P;P h ...,P n ) 



lim 

<3i ^-Pj, 



/or distinct Qi's (the order in which the limits are taken is irrelevant). 

Proposition 4.4. For given Pi, . . . ,P n e X g , we find Qi, . . . ,Qi with i = g for n > g 
and I = n otherwise, such that 

Pi + P 2 + --- + P n + Q 1 + Q 2 + --- + Qe-(n + l)oo~0 

by taking the zero-divisor of f-i n (P] Pi, ■ ■ ■ , Pn)- For each Qi = (xi, yi), by letting —Qi = 
(xi, —y^, we have the addition property, 

Pi + P 2 + --- + P n -noo~ (-Qi) + (-Q 2 ) + -.. + (-Q e ) - (-£)oo. 

As usual, the symbol of addition in the (free abelian) divisor group is used as well for 
addition up to linear equivalence. 

Remark 4.5. We should note that the hyperelliptic involution 1 : [x, y) (x, —y) induces 
the [— l]-action on J g , defined by u 1— > —u. 
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4.3. The analytic addition formula. We have the following addition formula for the 
hyperelliptic o functions (Theorem 5.1 in [?]): 

Theorem 4.6. Assume that (jn,n) is a pair of positive integers. Let (xi,yi) (i = 
1, . . . , m), (Xj,y'j) (j = 1, . . . , n) in X g and u G K~ 1 (W m ), v G /c -1 (W n ) be points such 
that u = w((xi,yx), . . . , (x m ,y m )) and v = w((x[, y[), . . . , (x' n ,y' n )). Then the following 
relation holds: 

nl=0 *m+n((xi, yi), . . . , (X m , y m ), (x[, (-lYyi), • • • , «, (-l)Vn)) 



(4.2) =5{g,m,n) 



x 



1 



, ._ 1 f^2{{x h y i ),{x' p y' j )) 
where 5{g,m,n) = (_i)sn+in(n-i)+mn_ 

Theorem 14.61 with m = g and n = 2 leads to the following Corollary: 
Corollary 4.7. Let (x i ,y l ) G X 5 (i = l,...,g), (x / j ,y' j ) G X 9 (j = 1,2), u G C g , 



v :— v 



|L1 + v [2] G k _1 (W2), and v [j] G kT^Wi) (j = 1,2) be points such that 



^((^1; Z/i)> ■ ■ ■ j (^gj and = w((x'j,y'j)), (j = 1,2). T/ten the following relation 
holds : 



(4.3) 

where H(it, u) zs e^ita/ to 



<r(it + f)cr(-u — f) 
cr(it) 2 (Jb(v) 2 



^ (xj x^)(a;j X2)F'(xi) 



anct! F(x) := (x - xi)(x - x 2 ) • • • (x - x s ) = li g ({x, y)\ (x u yi), . . . , (x s , y g )) and F'(x) 
dF(x)/dx. 

Proof. By letting A(x 1; x 2 , . . . , xi) be the Vandermonde determinant, i.e., 



(4.4) 



A(xi,x 2 , ...,x e ) 



1 X\ 

1 x 2 



X i 



.X o 



X o 



h3=hi<3 
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we have 



* fl+2 ((^i, Vi), {x g , y g ), K, ±y[), (4, ±y' 2 )) 



(4.5) 



1 X\ 

1 x 2 



1 Xg 

1 z'i 
1 x' 2 



x\ y 1 
x 2 y 2 



J9 _!_„,/ 



x\ ±y l 
i 

2 



x'% ±y' 



, . . . , J/j, . . . , Jjg, Jj 



i=l 



i=l 



Thus for the case m = g and n = 2, the right-hand side of the formula in Theorem 14.61 is 
equal to □ 

Baker |B1| §11.217], \B3\ p. 138] proves the following: 

Proposition 4.8. Let (xi, yi) G X g {i = 1, . . . , g) andu G C 9 such thatu = w((x\, yi), . . . , ( 
The following relation holds for generic x\ (i = 1, 2), 



(4.6) i=1 i =1 



where 



{x^ Xi){x 2 Xi)F' (^Xi) 



f(x' 2 )F(x[: 



+ 



f(x 1 , x 2 



(x[ - x' 2 ) 2 F(xi) (xi - x' 2 fF{x' 2 ) (x[ - x' 2 ) 



I \2 ' 



/(xi,x 2 ) = E x i x 2(-Wi(a:i +x 2 ) + 2A 2i ). 

Corollary 14.71 and Proposition 14.81 yield the key proposition in this article. 
Proposition 4.9. For the variables in Corollary \4-7\ the following relation holds: 



(4.7) 



a(u + v)a(u - v) _ f(x[, x' 2 ) - 2y[y' 2 _ ^ ^ 



a(u) 2 a\>(v) 



(xi - x' 2 ) 2 



, / i— 1 m— 1 



i=l j=l 



This corresponds to Fay's formula [Fj (39)], which is the basis of "Fay's trisecant iden- 
tity". 

Remark 4.10. For the g = 1 case, t^ll ^ — .^f 1 ^ 2 = p(^ 1 + v 2 )- Thus if we let v 3 = 

(xi-x 2 ) 2 

V\ + f 2 and X3 = p{vs), (14. 71) recovers (12.41) . 
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Corollary 4.11. For the variables in Corollary \4-7\ with t>M = v^ 2 \ we have 

2«"')<t( m -2i.I") 3 " 

(4.8) *{u)W*M)' tf}3 



0~\U + ZV l ' )0~(U — ZV l ' I , Y^V^ ( \ li+j- 



HI, 



— lim E(u,v). 



where 



h?{x) := |^ - fl 2 (x), fUx) := ^A^x 2 *" 1 + i(i - 

Belokolos, Enolskii, and Salerno gave the following relation [BESj Theorem 3.2]: 

Corollary 4.12. For the variables in Corollary \4-7l with = or (x^y'z) = oo, we 
have 



v(u + v [1] Mu - yW) _ /g _^ , 



=i 



= F(x' 1 ) = (x' x — x%)(x' 1 — x 2 ) ■ ■ ■ (x[ — x g ) 
= Hg((B f uVi)\ ( x i, yi), • • • > (xg, y g )). 

Proof. We divide both sides of (14 .7p by x' 2 ~ x . By taking the appropriate limit, we obtain 
the equality. □ 

Remark 4.13. We conclude this Section by elaborating on the geometric properties 
of the a function which make the addition formulas particularly suited for integrating 
equations of dynamics. As stated in the Introduction, originally the sigma function was 
discovered by Weierstrass [W] in order to express a symmetric function of the points of 
a hyperelliptic curve, which he called "al" function, in terms of rational functions. The 
al function is a root function, equal to a/ F(b a ) for a zero b a of f(x), giving one of the 
branch points of X g . The focus was on constructing a dictionary between abelian and 
rational functions. In particular, a goal was the solution of the "Jacobi inversion" , namely 
returning the symmetric functions of the divisor from a point on the universal cover of 
the Jacobian, as in the genus-2 case: p 22 = x\ + x 2 , <pi\ = x\x 2) where the points of the 
divisor P\ + P 2 are Pi = (xi,yi), i = 1,2. (see also Remark l5TTT 2)(b)). The most natural 
application is then the explicit realization of the group structure of the Jacobian (the 
generalized Frobenius-Stickelberger relation in Proposition I4.2[ which gives the addition 
structure, shows a simple connection between the affine coordinate ring R g and the group 
law on the Jacobian J'g), reflecting addition in the free abelian divisor group, and this is 
achieved in Corollary I4.11[ vis-a-vis Proposition 14.41 
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5. The cx-function solution of the Toda lattice 



In this Section, we give the a-function solution of the Toda lattice equation. 

First, we identify some algebraic identities that hold for vector fields on the Jacobian. 
Vector fields on the Jacobian are understood to be translation invariant; equivalently, 
they are elements of the tangent space at the origin. It is important that we use algebraic 
functions on the curve to write their coordinates in the abelian variables (ui), but in 
doing this, we make the convention that we are on a suitable coordinate patch on the 
Jacobian, where the Abel map from S 9 X g can be inverted; as explained in |Mu[ §3], there 
are choices involved and one has to check that the vector field in question is well defined. 
Our formulas would hold on this suitable affine patch anyway, since outside it, the Toda 
orbits become of smaller dimension, as mentioned in the Introduction; that case will be 
addressed in our forthcoming work. Moreover, as usual, we view the vector fields as 
derivations on the universal cover, namely on the ring of functions in T(C 9 , 0(C 9 )). 

In order to give our solution of the Toda lattice equation, we need only two vector fields 
on the Jacobian, each associated in the same way to a fixed point (x'j, y'j) G X g , j — 1, 2; 
however, in order to relate our solution to Baker's differential equations (Remark 15. 7p . 
we distinguish two further fixed points (x^y^) G X g , j = 3,4, for which we use Baker's 
expression of the derivatives in terms of algebraic functions on the curve as opposed to 
abelian coordinates on the Jacobian; the dictionary between the two expressions is given 
in Lemma 15.21 (b) . 

Definition 5.1. For (x h y { ) G X g (i = 1, . . . , g), u = w^x^y-i), (x g , y g )) G C 9 , and 
(x'^y'j) G X g (j = 1,...,4), we let 




Lemma 5.2. Let (xi,yi) G X g (i = l,...,g), (x'j, y'j) G X g (j = 1,...,4), u G C 9 , 
and G /c _1 (VVi) (j = 1,...,4) be points such that u = w((xi, yi), . . . , (x g , y g )) and 
i>b'l = w((x'j, y'j)), {j — 1, • • • , 4). We have the following expressions: 
(a) For each j = 1, 2, 3, or A, 

2yid xi 
2y2d X2 



o 



Dj A(xi,x 2 ,...,Xc 



5-/ TPI 



yiF(x', 



x L 
x 2 



d 
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9-1 

9-1 
x 2 



X 



and for each j, f — 1, . . . , 4, 

[Dji Dj>] = DjDji — DjiDj = 0. 
(b) For vv> = w(x'j, y'j) with j — 1, . . . , 4, 

9 1 <A 9 



(c) For h G r(C , 0(C 9 )) and j = 1, . . . , 4, 

+ v w ) := 2y' j ^- r h(u + v {j) ) = D 3 h(u + v {l) ). 

Proof. A direct calculation gives the results. □ 

Lemma 5.3. For the variables in Corollary \4-7\ we have 
(a) 

D x log a(u + v) = - (D 1 log E(u, v) + Dy log E(u, v)) + D 1 log a(u) + Dy log cr^u), 
(b) 

DiD 2 log cr(u + u) = -D\Di log 5(«, f ) + -DiD 2 ' log t>) + D X D 2 log a (it) 

= -DiDalogS^u) + ^£> 2 £>v logS(u, u) + Di^log^u). 

Proof. By taking the derivatives of the logarithm of (14. 3p . we obtain 

Di log H(iz, f ) = Di log cr(w + f ) — 2Di logcr(it) + loga(u — t>), 

Dy log H(it, u) = Dy log cr(w + v) — 2/}]/ log cr|,(f ) + Dy log ct(m — f ), 

and Lemma 15.21 (c) yields the claims. □ 

Now we can give the sigma-function solution of the Toda lattice equation: 

Theorem 5.4. Let (a^jfc) G X g (i = l,...,g), «,^) 6l 9 ue O, and «W G k _1 (>Vi) 
6e points such that u = w((xx, y±), . . . , (x g , y g )) and t>M = w{{x' x ^y'^}). We define c := 
2vM A = <7 b (c)L> 1; 

V(«) := V(«, «W) := E PyW^" 2 ' V ^ c ) : = /i,2K)» 

i=l i=i 

and i := (t n , t 12 , . . . , G C 9 

i=\ J 00 

TTien with the coordinate change u = nc + t 1 - + t (which defines t x ), we have 
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(1) 

-Dl\og{V(t + nc + t^)-V c (c)) 

= V(t + (n + l)c + t ± )-2V(t + nc + t ± ) + V(t + (n-l)c + t ± ). 

(2) The Hirota bilinear equation, 

a(t + nc + t L )Dla(t + nc + t 1 ) - D x a{t + nc + t L )Dia{t + nc + t 1 ) 
{) " ] -V c (c)a{t + nc + t^f - ait + (n + l)c + t L )a{t + (n - l)c + t x ) = 0. 

(3) The Toda-lattice equation, by letting V n (t + t ) := — V(t + nc + t -1 ) and q n (t) := 
-log(V„(t + t ± )-V c (c)) ; 

(5.3) -Dlq n (t) = e- qn+1 - 2e~ qn + e -9 - 1 . 

Proof. The claims follow from Corollary 14. 11[ by direct verification, as in the genus-1 case, 
using the definition of the vector fields. □ 

We also have the sigma function solutions of the two- (time-) dimensional Toda lattice 
equation JEH 3.5.1, (3.107-108)], d 2 q n /dt 1 dt 2 = V n+1 - 2V n + K-i, q n = log(l + V n ), n 
any integer: 

Theorem 5.5. Let {x hVi ) E X g (i = l,...,g), (x^yfi E X g (j = 1,2), u E C», 
v ■= t>W + t>M 



G K 1 (W2), ond E K 1 (Wi) (j = 1,2) 6e points such that u 



w{{xx,yx), . . . , (x g ,y g )) and = w((x' p y',)), (j = 1,2). We define c := v [1] + v [2] , 



9 9 



/ i— 1 / j— 1 



i=l J = l 



V C („M „ M ) ;= 2^-/K,4) 



37 1 3/ n 



and i,- := (t,i, tj 2 ,..., tj g ) E C 9 with 

9 r(%i,yi) 

*ifc:=(4) 1_fc 2 / 0' = 1 > 2 , and* = 1,2,..., y). 

i=i ^ 00 

T/ien with u = nc + t 1 - + ^ + t 2; we /iai>e 

-DiD 2 log(V(nc + h + t 2 + t x ) - V c (c)) 

= V(ti + t 2 + (n + l)c + t x ) - 2V(ti + t 2 + nc + ^) + V(*i + t 2 + (n - l)c + ^). 

Proof. The claims follow from Proposition 14.91 by direct verification. □ 

Remark 5.6. As Hirota notes, a two- (space-) dimensional version of the Toda system 
ought to have two independent discrete variables, but at the time, he had only found 
solitons for the two-times case [Hi| §3.5.1, Remark]. Instead of Theorem I5.5[ as a different 
type of generalization of Theorem 15.41 involving the sigma function on the Jacobian, 
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implementing a second spatial (discrete) variable would entail choosing another G Wi 
say, and adding to the argument of the cr-f unction a vector mc\ +nc2, with Ci = 2t>M and 
C2 = 2^W, and "considering a two-dimensional version of the force term on the right-hand 
side" of (1) in Theorem 15.41 [ibid.]. We would need a different kind of addition formula, 
which we have not yet developed. 

Remark 5.7. We would like to stress that a single formula underlies the validity of the 
integrable hierarchies: the differential operators that Baker used to obtain what we call the 
KdV and KP equations are always of the type Di we defined, namely a linear combination 
of the Abelian coordinates given by the basis of holomorphic differentials which involve 
the equation of the plane curve, against coefficients that are powers of the x-coordinate of 
one point of the curve; generically, distinct points yield independent vector fields. Specific 
to our situation: 

(1) Corollary Baker's formula glj) in [EQ p. 328] and [B3l p. 138], and Fay's 
formula (14. 7p in |F, (39) in p. 26] are essentially the same. We notice the following 
two facts: 

(a) Baker derived the KdV hierarchy and KP equation by using his formula (14. 6p . 
using the vector fields D 3 and _D 4 , cf. Definition 15.11 |Mat2j . 

(b) Fay derived his famous trisecant identity, which is equivalent to the KP hierar- 
chy, by using formula (14. 7p . which is the hyperelliptic sigma-function version. 

(2) The following relationships hold among the formulae: 

(a) The two-dimensional Toda equation in Theorem 15.51 is the same as formula 

IZD. 

(b) When = in Theorem 15. 5\ we obtain the Toda lattice equation of Theo- 
rem |5]H and when = 0, the Toda equations are obtained by differentiating 
F(x[) (Corollary I4.12p . The function F(x[) is a polynomial of degree g in the 
variable x[. By applying the vector field D 3 , we obtain the "Mumford triple" 
(which is called (U, V, W) in |Muj . three polynomials that parametrize the Ja- 
cobian outside a theta divisor (excluded are certain g-tuples of points (x^, yi) 
in special position). We note that the coefficients of Mumford's (U, V, W) 
involve only the abelian functions p g i, i > [Mu, §10]. The KdV hierarchy 
follows again using formula (14. 7p . 

(3) We stress again the connection between algebraic and abelian functions: the dif- 
ferential operator Dj has an expression that involves only (and acts upon) the x^s 
(cf. Lemma [5.2p . On the other hand, (14. 3 p is also given by the affine coordinates 
of points of X g . Hence the Toda lattice equation is an identity defined over S 9 X g . 
It can also be regarded as a relation among the functions over J g and S 9 X g . 

Remark 5.8. We conclude by comparing our solution to some of the methods that were 
used to obtain algebro-geometric Toda flows. 

(1) van Moerbeke |vMj . following his work with Kac on the Toda lattice and Jacobi 
matrices, reported below in Section El gave a description of the isospectral flows in 

16 



terms of linear flows on a Jacobian, and with Mumford |vMMj . the algebraic coor- 
dinates for the flows and the solutions in terms of theta functions. The flows are 
linear combinations of, essentially, the defined above, for (x'^y'j) a branchpoint 
of the hyperelliptic involution. Fay's trisecant identity and its derivatives along 
the Di are used |vMMt §5, Proposition] to show that the flows are Hamiltonian 
vector fields that preserve the spectrum of the matrices. 
(2) Algebro-geometric solutions to the Toda lattice can be found in [GHMTJ. These 
authors in their extensive work also used the spectral curve of the tridiagonal 
matrices whose deformations, equivalent to the Toda lattice, are given below in 
Section O To solve them, using the divisors given by auxiliary spectra and via 
eigenvectors expressed in terms of theta functions, the authors work out a "discrete 
Floquet theory" (as had done Kac and van Moerbeke) and solve the "Dubrovin 
equations" on the expansion of the Green's function by Abelian functions. 



6. Periodicity of the Toda-lattice solution 

In this section we turn to the problem of periodic solutions of the Toda lattice. Spatial 
periodicity is of physical interest, in the lattice case, so the requirement amounts to finding 
a point of finite order N on the hyperelliptic curve: c = 2w((x[, y[)) such that for u € C 9 , 

u + Nc = u modulo A. 

Hyperelliptic curves that admit such a point are special, and were called "Toda curves" by 
McKean and van Moerbeke |McKvM] . who proved that they are dense in moduli space. 
For the same reason, points of finite order in the Jacobian that come from the curve, 
or from the sum of a specific number of points on the curve, give periodic orbits in the 
billiard (completely integrable Hamiltonian) system in the ellipsoid |DRj , and are related 
to Poncelet's closure (cf. Appendix). 

The finite-point condition was investigated by Cantor and Onishi [oj E] using the 
division polynomial ip2N, an element of R g . Similarly, we investigate the periodic solutions 
of the Toda lattice. 

6.1. Division polynomials ip2N- Traditionally, division polynomials F n (x,y) arise in 
expressing the coordinates of nP in terms of those of P, a point of an elliptic curve in 
Weierstrass form. In particular, given our present interest, we call "the division polyno- 
mial" an element of the ring of functions on the affine curve, whose solutions are points of 
finite order in the Jacobian (where the Abel map is, as usual, based at oo). The division 
polynomial for the genus-one case was studied by Kiepert |Ki] and Brioschi [Br] . We 
call Kiepert-type polynomial and Brioschi-type polynomial the genus-one version of the 
division polynomial ip n , a polynomial whose zeros P satisfy the condition nP = 0, more 
precisely, nw(P) = modulo A. 
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Referring to |6| Definition 9.2], the hyperelliptic version of the ip n function for genus 
X g over w(X g ) = k _1 Wi is defined by 



a^n {nu) 



A zero 0) of tp n has the property that nu G K~ l W g -i. The transformation law 
under translation (13. 3p allows one to check that ip n belongs to R g . Thus it is a natural 
generalization of the classical Kiepert formula, or the division polynomial. By taking 
limits of Proposition 14.21 along the Abelian variables, we can give an expression for if) n in 
terms of 0j's in R g [oj Theorem 9.3]. In |MP1] . we proved the following: 

Theorem 6.1. Let n > 1 be a positive integer. For 



"n,g 



l!2! • • • (n - V 



C01 


002 


Otpn~l 


dui 


dui 


dui 


d 2 ^ 


d 2 02 


d 2 <Pn-l 


du± 2 


dui 2 


dui 2 


a™- 1 01 


d n 02 




dux"- 1 




&U!™- 1 



with ipi = 1, the vanishing of ip n at a point P of the hyperelliptic curve X g is a necessary 
and sufficient condition for uj(n- P) to belong to W fl _i, where u is the Abel map uj : X g — > 
J g . Here e n ^ g is a plus or minus sign. 

Further letn(> g), k(< g) and I := g—k — 1 be non-negative integers. For a hyperelliptic 
curve X g , the vanishing of ip n+l , . . ., ip n+ i, ip n! ip n -i, ■ ■ ■, ip n -i> a t a point P of X g is a 
necessary and sufficient condition for u(n ■ P) to belong to V\4. 



Cantor gave a Brioschi-type expression for the ^-function [0, IMatl] , 

e' n j2yr^W ■ T^f_ 1)/2 (y, £) for n £ g mod 2 

4, 3 (%) n(n - 1)/2 ■ T(n-lUy> s) for 71 = 9 mod 2 - 

Here e' is a plus or minus sign and T n m ^ is a Toeplitz determinant |Matl] , 



Ti m) I g(8) 



d_ 

ds 



9 



m+n— 1] 



9 



m+n] 



[m+n— 2] 
[m+n— 1] 



g[m+2n-3] g[m+2n-4] 
g[m+2n-2] g[m+2n-3] 



m+1] 
m+2] 



[m+n-1] _ 



[m+n] 



9 



g[m+l] 

Vn+n-2] 
m+n— 1] 



d 



and Ti m) 

g(s) is a function of an argument s and 



g(s), — 1=1 when m is a non-negative integer and £ is a negative integer, 
ds 



1 d n 

n\ ds' 1 
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■9{s). 



We note for y 2 = f(x) that y 2n l d n y/dx n is a polynomial in x and coprime (in the 
sense of not vanishing together on a point of the curve) to f(x) in general. Hence the 
function y^ 2 " 1 ^- 3 ^^ (y , f- ), that is, 



dx ' 

y2m+2n—3 y[m+n—l] y2m+2n— 5y[m+n— 2] 
y2m+2n—ly[m+n] y2m+2n—3y[m+n—l] 



y2m+An-7 y[m+2n-3] y2m+An-9y[m+2n-A] 
y2m+4:n-3y[m+2n-2] ^2m+4n-7 y[m+2n-3] 



y2m+ly[m+l] y2m-ly[m] 
y2m+3y[m+2] y2m+ly[m+l] 



y2m+2n—3y[m+n—l] y2m+2n—5y[m J rn—2] 
y2m+2n+l y[m+n+l] y2m+2n—3y[m+n—l] 



is an element of C[x] and coprime to y 2 = f(x). Hence ip n (u) can be expressed by 

{(2 2/ )9(9+ 1 )/ 2 a n ( x ) for n - g = odd, n > g + 1 
{2y) 9 ^ 9 ~ x ^ 2 a n {x) for n — g = even, n > g + 1 
(2y) n ^ n ~ r> / 2 a n (x) otherwise, 

where ct n (x) is a polynomial of x and coprime of y for n > g + 1, and a n = for 
n < # + 1. As shown in [C], the degree of a n (x), {n > g + 2) is 

g(n + g)(n- g) - g(2g + 1) 



2 

g(n + g)(n-g) 



deg(a n 

2 

Definition 6.2. VFe define 



for n — g = odd 
for n — g = even. 



$ n := {P G X 9 | P is a zero o/a n }, 

an<i /or n > g, 

E n := $ n - g +i n . . . n $ n -i n $ n n $ n+i n . . . n 

One should note here that there is no guarantee that S n is not empty. 

Theorem 6.3. Let 2N > 2g + 1. For a hyperelliptic curve of genus g which has a point 
(x'i,y[) G H 2 at, V(w) m Theorem\5.J\ is the periodic solution of the Toda lattice equation 



such that V{u) = V{u + Nc) with c = 2w(x[, y[). 

Remark 6.4. (1) In the g = 1 case, since an elliptic curve is a divisible group, there 
exist a point {x'^y'^) which is a zero of a ip^; as mentioned in Remark I4.10[ for 
every point (x'i,y'i) we have a point (x^y^) = 2(x r l ,y[). Thus, V in Theorem 15 .4} 
is a periodic solution of the Toda lattice equation, provided V{u) = V{u + Nc) 
with c = u>(4,2/2)- 

(2) We do not address in the present work the important problem of finding real- valued 
solutions V{u). 
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Example 6.5. Case g = 1, JV = 3 and N = 4. 

We consider the elliptic curve y 2 = x 3 — x (which has an extra automorphism of order 
two - a property that is not related to points of finite order but is usually accompanied by 
a large number of exact solutions for the coordinates of points on the curve that satisfy 
algebraic conditions). The division polynomials are given by 

^i = l, 
^2 = -2y, 

?/> 3 = My 2 (x A -6x 2 - 1), 
^ 4 = ~2y{x 2 + l)(a; 2 + 2x- l){x 2 -2x- 1), 
V> 5 = 32x 14 - 187a; 12 - 64a; 11 + 2a; 10 + 320a; 9 - 233a; 8 
+ 320a; 7 - 52a; 6 - 64x 5 - 61a; 4 + 50a; 2 + 1. 

For x'n = (1/3)^9 + 6^3, we have N = 3 and for x' 3 = a/2 + 1 or a zero of ^4 we have 
N = 4. 

Lemma 6.6. Let 2N > 2g + 1. If P := (x^y^) G has the property that IP are 
distinct for different £ G {1, 2, • • • , 2iV} ; then 

(X - x[)(x - x' 2 ) ■ ■ ■ (X - x' 2N )\lp2N+m(x, y) 

where (x' e , y' e ) := &(x , 1 ,y' 1 ) and m = — g + 1, . . . , 0, . . . , g — 1. 

Proof. By assumption, ±£P are exactly the points of the set S 2 at. □ 

6.2. A periodic Toda lattice. In this subsection, we consider the relation between 
Theorem 16 . 3 1 and an algebraically completely integrable system^ originally studied by Kac 
and van Moerbeke |KvM] . Using the solution given in Section [51 we find the explicit 
form of Flaschka's coordinates for the Toda lattice. The Hamiltonian of the Toda lattice 
equation is 

. N N 
k=l k=l 

where Pk = Pu+n and Qk = Qk+N- For Flaschka's coordinates, = exp(Qk — Qk+i) and 
bk = —Pk, the equations of motion under H become 

— 0^ = ak(bk+i — bk), 

(k = l,2,...,N). 

d h 

T7°fc — a k — a k-l- 

at 

Remark 6.7. The Toda Hamiltonian system admits the time inversion 1 1— > —t, and this 
corresponds to the hyperelliptic involution on X g . 



5 Different systems have been variously referred to in the literature as "periodic (generalized) Toda 
systems" ; for extensive information on definitions and solutions we refer to the monographs [AvMVl IV] . 
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For brevity, we introduce the notation 

<7(»)(i; t x ) : = a (t + nc + t 1 ), a (c) := a b (c), 

C {n \t;t L ) := ^ a ,r 1 C i (t + nc + t ± ), C (C) := \d v log ^ (c) , 



i=i 
s 



^(t;^) := £ x' 1 ^'- 2 ^(t + ^c + t ± ), := / li2 (4). 

Proposition 6.8. Using the sigma-function solution of the Toda lattice equation in The- 
5.4\ (2), the Flaschka coordinates for the Toda lattice are expressed as follows: 



orem 



V-i = A log Jn-ifff) - Cc = C w (t; t x ) - (^(t; t 1 ) - C (c) . 

Proof. By definition of Pfc = — bk, 

Pn-i ~Pn = & +1) (t; t 1 ) - 2C (n) (t; t X ) + C (n-1) (i; 



p 2 - p 3 = c (4) (*;^) - 2C (3) (t;t ± ) + C (2) (^;^), 
p - p 2 = c (3) (t; ^) - 2C (2) (t; ^) + C (1) (*; 

Po - P = C (2) (t; t^ - 2C (1) (t; t x ) + C (0) (^; 



with 

-P n = ( {n+1) (t; t^ - C {n \t; t*-) - (C (1) (t; t X ) - C { °\t; **-)) - P . 

Since the total momentum should be invariant, we set 

P = -(C (1) (t;t ± )-e\t;t ± ))+ Po . 

where p$ is a constant corresponding to Q c \ Then the equation db^/dt = — a^-i holds 
by Lemma [5.31 and Corollary 14. 7[ equation da^/dt = ak(bk+i — bk) by Theorem 15 .41 □ 

Proposition 6.9. The a n 's and b n 's given in Proposition ^.^ are rational functions of Xi, 
V% (* = 1, • • • ,9) a nd x' x , y[. 

Proof. The translation law (13.31) for the functions shows that a n and b n are meromor- 
phic functions over the Jacobian J g ] Lemma [5.21 gives b n as a meromorphic function over 
S 9 X g x X g , derived using algebraic vector fields whose coordinates are rational functions 
of x i; yi (i = l,...,g), x[, y[. 
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On the other hand, using Theorem 14. (j\ a n is given by 

a(u + 2nv [1] + 2v [1 ^)a(u + 2nv^ 1] - 2v [1] ) 
<x(u + 2nvW) 2 a ir (2vW) 2 



lim 



(* s+2 „((Xi, J/i) ... , (4„+2) ?/2n+2))^2(K, J/£)> 04, 2/ 2 ))) 



2 



9 2 2n+2 2 



Here ((xi,yi), . . . , (^2n+2) 2/20+2)) m both numerator and denominator of the equation is 
an abbreviation for ((xi, yi), . . . , (x g , y g ), (x' 3 , y' 3 ), ... , (x' 2n+2 , 2/20+2))- This completes the 
proof of the statement. □ 

Remark 6.10. To exemplify Proposition 16.91 we write b n in the g = 1 case: 

C{u + v) -((u)-((v) 



x(u) — x(y) 

Now we give an inverse of Theorem 16. 3} together, they constitute our main result 
relating the periodic Toda system to the sigma-function solution: 

Theorem 6.11. Let2N > 2g + l. If a hyperelliptic curve of genus g has a point (x^y^) G 
X g \oo such that V(u) in Theorem \5.J\ is the periodic solution of the Toda lattice equation, 
i.e., V(u) = V(u + Nc) with c = 2w(x[,y' 1 ), then {x'^y'^j belongs to S 2 at. 

Proof. The periodicity condition V(u) = V(u + Nc) is equivalent to the periodicity of a n 
and b n due to Kac and van Moerbeke |KvM] . The expression for the Flaschka coordinates 
in Proposition 16.81 allows for a prefactor e P u of o~(u), (3 a constant vector in C 9 , in terms of 
which the periodicity of a n and b n , (a n = a n+ N, b n = b n+ N), is equivalent to the equality: 

(6.1) e 2Nt(3v[1] a(u + 2Nv ll] ) = a(u) 

for every ±u G Wk (k = 0, . . . ,g), where = w(x' 1 ,y' 1 ) ^ (by assumption). Noting 
2N — (g — 1) > g, this implies in particular, setting u = ±kv M ^ 0, 

e 2iv t ^W (7 ((2iV ± k)vW) = a(±kvW). 

The right-hand side vanishes for (k = 0, ...,g — 1) since u = ±/ct>M belongs to the 
(translated) theta divisor, whereas o~$(v M) does not vanish. 

Moreover, we recall that a point (x,y) in X g with y = is such that 2w(x, y) G W g -i, 
which still satisfies the conclusion since N > 1. Hence (16.1 j) implies that (a^,?^) G 

Remark 6.12. Given the connection between Theorem 16.31 and theory of Kac and van 
Moerbeke provided by Theorem I6.11[ we note in addition: since it is known that for 
certain multi-indices 7 = (ji, ...,j g ) and for ±u G re -1 VVfc (k = l,...,g — 1), letting 
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<9 7 := d^l ■ ■ ■ dZ 3 g , the derivative <9 7 cr(u) vanishes [0 [ IMP2] . by differentiating the equation, 
the a function satisfies (k = 1, ...,£) for a suitable £, 



d^a^N - £)v [1] ) = <9 7 cr((2N - £ - 



<9 7 cr((2N + £)v 



0. 



As we have assumed that i>M G does not vanish, the vanishing of a and its derivatives 
on multiples of v [1] is a condition of flag- variety type (cf. |A"M| IAHM| EHJ EH iKo] ). 
Studying the topology of these Toda orbits of smaller (than generic) dimension was the 
original motivation for our work, which we plan to use for a sequel to this paper. 

6.3. Hyperelliptic curve X g ^^\ for the periodic Toda lattice. The Lax matrix for 
the periodic Toda lattice is now given by 



(b x 1 

a x b 2 

•■■ 

\w ■ ■ ■ 




1 

a N-2 &JV-1 
a N-l 



a N w 1 \ 




1 

b N 



J 



The characteristic equation for £ defines the hyperelliptic spectral curve: 

nN 
i=l° 



det(£ - z) 



w 



w 



V(z) 



0. 



which gives the affine curve (w, z) of genus N — 1, 

N 

(6.2) ^ g ,iv-i : w 2 -V{z)w + Y[di = 0. 

Here P is given by 



i=i 



where 



V{z) :=A^ 1>N (z)-A^ 2!N ^{z), 



bin % 



A 



per 



d m ^m+1 Z 1 













a n~2 b n _i — Z 1 

a n _i b n — z 



V{z) := (-l) N z N + N(-l) N+k l k z 

k=l 
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The coefficients of the powers of z are given by 

N N 
i=l i>j i=l 

(6-3) 

N 

X N = A pcr li7V (0) - a iv AP er 2jJ v_ 1 (0), X N+1 = J]^- 

i=i 

We refer to X 9t N-i as the periodic Toda curve. Since the characteristic polynomial det(£- 
z) is invariant under the Toda flow, these coefficients give the Hamiltonians, i.e. 



(6.4) ]ft Ij = ° U = 1 i 2 t--iN + 1) 



The set {Xj : j = 1, . . . , N + 1} gives an involutive, complete family of integrals of motion, 
which guarantees the complete integrability of the Toda lattice. In particular, X\ and 
Xat + i can be expressed as follows: 



Example 6.13. From the formulae of (a n , b n ) in Proposition 16.81 the following expression 
for two Hamiltonians in terms of Abelian functions follows directly: 

N 

Zi = J> = N( (c \ X N+l = JJo, = (a^)- 2N - 

i i=l 

Lemma 6.14. Xj (j = 1, .. .,N+ 1) can be expressed as rational functions of 
only. 

Proof. Condition (16.41) shows that Xj (j = 1, . . . , N + 1) does not depend upon (xj, yj), 
j = l,...,g. ' □ 

Proposition 6.15. The hyperelliptic curve X g ^-i of §6. 2$ admits a morphism to the 
curve: 

2N N 

w 2 = l[(z-z l )=V(z) 2 -4l[a i , 

i=l i=l 

with w := 2w — V(z). The coordinates of the Weierstrass points of X g ^-i are rational 
functions of x[ and y[ . 

Proof. The Weierstrass points have ^-coordinates which correspond to w = 0. The co- 
efficients of the corresponding polynomial in z are rational functions of x[ and y[ by 
Lemma I6.14[ and the fundamental theorem of elementary symmetric functions gives the 
conclusion. □ 
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Example 6.16. For g — 1, we consider again the curve y 2 = x 3 — x or y 2 = 4(x 3 — x). 
Then, 

(3x4 2 - x - x' 3 - 2yy 3 - x' 3 ) 



ai(i) 



(x-x 3 ) 2 



a 2 (t) = 2(x — x' 3 ) 



(6.5) , . o . , . / (y' 3 4 + 8^' 3 3 + (-36xx' 3 2 + (2(12 + 9x 2 ))x' 3 - 12x)y' 3 



-3x' 3 3 + x' 3 + 2y' 3 - 2yy' 3 + 3xx' 3 2 



x 



T 2 _|_ q T ' 6 _ fir' 4 -I- r' 2 -I- 1 9tt' 3 — 9rr' — 1 Srr' 5 — fir 2 r' 2 -I- Qt 2 t' 4 



34 3 + x 3 + 2y 3 2 -2# 3 + 3x4 2 



X 



(6.6) 6o(t) = ^-%, 6i(t) = ^ %, fc(t) = 



/ J J- V / / ' z V / / ' 

X — X 3 </a;,l — ^3 7a;,2 ~~ x 3 

where (773,1, %i) and (rj Xt2 , ^,2) are the solutions of the equations /ii((<2, w); (x, y), (x 3 , y' 3 )) 
and fjLz((ZjW); (x,y),2(x 3 ,y 3 )) = in the variables (2,10), with: 

// w -w / ./u (x> - y 3 z + z;^ - xw + x£ 3 - ;yx 3 ) 

Mi((^w); (^3/), (^3^3)) : = ^7—^ > 

l x 3 X ) 

si *s ~i 1 -/xx —zx 2 + xz 2 + 2xx'^wy^ — 2zyx',,y'n—zx'J 2 — 3xz 2 x'J 2 
fjL 2 ((z, w); (x, y), 2(4, y' 3 )) := 3 ,f f3 y 393 3 3 - 



x - 3xx 3 + x' 3 + 2yy 3 + x 3 
^ 3zx 2 x' 2 + xx' 2 + ra', 4 — zx', 4 — 2wx 2 v'^ -\ 



3zx 2 x' 3 + xx' 2 + iXj 4 — zx' 3 — 2wx 2 y' 3 + 2yz 2 y' 3 + z 2 x 3 3 — x 2 x 3 3 + x' 3 z 2 — x 3 x 2 



x — 3xx' 3 2 + x' 3 + 2yy 3 + x 3 



6.4. Remark on a relation between X g ^-i and X g . Even though the Toda time flows 
have been linearized on both the Jacobians of the curves X g and X g ^_\ (cf. Remark 15.81) . 
and solutions have been expressed in terms of hyperelliptic abelian functions on each, the 
relationship between these curves is non-trivial. To name only one other example, the 
Kowalevski solution of the top has been compared by several authors with the Lax-pair 
solution, and the classical Arithmetic-geometric Mean has been shown to relate the (genus- 
2) curves |LM] . the Jacobian of the one is the quotient of the Jacobian of the other by a 
group of order 4. A different Lax pair provides a curve of genus 3, which covers an elliptic 
curve: the Prym variety of the cover is again isogenous to the genus-2 Jacobian [Mar] . 

In our case, we make some observations. The coefficients of A 5i at_i are given by hyperel- 
liptic (- and a- functions, which are transcendental functions of their arguments (points of 
the Jacobian). However, Propositions 16.81 and 16.91 show that they are rational functions of 
Xj's and y^'s; moreover, these coefficients are Hamiltonians of motion, so the t-dependence 
of the a n and b n fixes the curve X gt N-i, associated to a point on the Jacobian of X g . 
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We note that the cyclic group CV acts on the sets of a's and fe's by: a n h-» a n+ i 
and b n i— >■ b n+ \, via the addition on 2P (the image of the point P in the Jacobian has 
order 2iV): 2£P ^ (2£ + 2)P and {2£ + 1)P ^ (2£ + 3)P, and that the curve X g>N -i 
is invariant under this action, although this does not guarantee that X Si at_i admits a 
Cn action. There could be a relationship of Kowalevski type between the two Jacobians 
(which have different dimensions in general), such as quotienting by the group Cn- We 
plan to investigate this relationship, starting with small-genus examples. 

For example, in genus 1, the relation between X\ and X^jv-i could give another solution 
to Poncelet's closure problem (for Cayley's solution, cf. |DRt IGH] ): as we review in the 
Appendix, the periodic Toda flow corresponding to a point of order 2N on a given elliptic 
curve, which plays the role of our X g , also corresponds to a closed Poncelet polygon with 
2N sides; in this case, as we saw in subsection 16.31 the periodic Toda curve has genus 
N — 1, and could be viewed as an algebraic solution to the porism. 

We conclude the Appendix by a reference to another classical problem solved in terms 
of transcendental functions over an algebraic curve, which therefore could play the role of 
X g , a construction similar to ours yielding algebraic solutions over X 9t N-i- 

A. Appendix: Toda lattice and Poncelet's closure problem 

Poncelet's porism (cf. |DRj IGH] ) can be stated as follows: 

Theorem A.l. (Poncelet) Let C and D be two smooth conies in the real affine plane, such 
that C includes D. For an integer N > 2, if there exists a closed N -polygon inscribed in 
C and circumscribed about D, for every point P in C there exists a polygon whose vertices 
are in C and includes P, and segments are tangent to D. 

More generally, in the complex projective plane, the existence of such an A-sided 
Poncelet polygon corresponds to a point of order N on the elliptic curve determined by 
the two conies, and a point in the incidence correspondence of points of C and tangents 
to D. This interpretation in terms of a transcendental problem is one of the deeper ways 
to prove the theorem (cf. [G]), whereas Poncelet used elementary projective geometry, a 
subject that did not exist at the time. 

As can be expected, a point of finite order gives rise to many applications in the theory 
of periodic motion, and recently in |BZ] it was applied in a novel way to give a condition 
of Fritz John type on a Dirichlet problem for a planar domain bounded by an ellipse; in 
the same paper, the authors give the following explicit parametrization of the conies, and 
show that the vertices of a Poncelet iV-gon give a solution to the iV-periodic Toda chain, 
which we reproduce below. 

Without loss of generality, we assume the conic C to be given by the equation x 2 = yz 
and parametrize it by (x, x 2 , 1). Let the vertices of the Poncelet polygon be the iV-points 

{xf \ xf )2 , 1) (z = 1, . . . , N). Let D be defined by: 

(x,y,z)A t (x,y,z) = 0, 
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where 

/ ai a 2 a 3 
A = I 04 05 as 

Assume a 5 = 0. The dual conic D* of D is given by (X, Y, Z)A- n (X, Y, Z) = 0. A pair 
(P, L) e C x D* such that ?ei satisfies 

il + + zZ = 0, 

for P = (x, y, z) and L = (X, Y, Z). The relation is reduced to the elliptic curve E\ 

w 2 = — - — (x, x 2 , l)A\x,x 2 , 1), 
a 2 + a A 

where 

m= 7ab( Ma 4- Ml) 

and the h\, h 2 are polynomials in x. Poncelet's closure is equivalent to finding a matrix A 
as above and a point (x, w) belonging to E 1 such that it satisfies the equation of Kiepert 
and Brioschi, 

ip N ((x,w)) = 0, 

a criterion attributed to Cayley and proved in [GH] . We regard the equation of iPn((x, w)) - 
as the moduli equation for a given x. 

For such an A, Poncelet's theorem means that the vertex P n = (x n ,x^, 1) £ C {n — 
1, 2, . . . , N) satisfies the periodic Toda lattice, x n = p((n — l)u + t) \BZ\ §7.1], 

d 2 

(A 1) ~dP log ^ nu ° + t ">~ = ^ n + 1 ) u ° + *) _ 

- 2[p(nu + t)- p{u Q )\ + [p((n - l)w + t) - p(uo)], 

where 

(A.2) mo = 



/ (0) (0)n 

2w 



(in the previous notation, Xn ^ = p((n — 1)uq)). 

Lastly, we cite the problem of finding the general solution of the fifth-degree algebraic 
equation with Galois group PSL(2,F 5 ) = A 5 , alternating group on five elements. While 
the solution cannot be algebraic in the coefficients of the polynomial equation, classical 
authors such as Jacobi, Galois, and Klein |K1] . gave a solution in terms of the zeros of 
tps for an elliptic curve xE In the classical, fifth-degree case, Humbert |Hu] expressed 
the period-5 condition in terms of a Poncelet pentagon, but in addition proved that it is 
equivalent to the curve y 2 = {x—X\) • ■ ■ (x— £5) having real multiplication by the quadratic 



6 More recently, H. Umemura |Mui Chapter III.c] gave a solution of an algebraic equation of any degree 
n in terms of tlicta constants of a hyperelliptic curve of genus n, in other words, in terms of a Siegel 
modular function. 
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order of discriminant 5. Hashimoto and Sakai |HSj expressed the general condition for 
a hyperelliptic curve of genus 2, H 12 , y 2 = (x — X\) • ■ ■ (x — x 5 ) ■ (x — x 6 ) to have real 
multiplication of discriminant 5 again translating it into a condition for the existence of 
closed Poncelet pentagons. Mestre [Mej generalized the condition of real multiplication 
to genus g. Lemma [6.61 gives a criterion for the points of order 2N in X g in terms of the 
(genus-g) division polynomial; by analogy H\p might be replaced by a hyperelliptic curve 

Hg,N-X, 

V = (x — — x 2 ) ■ ■ • {pc — x 2 n)- 
Proposition 16.151 gives an algebraic relation between H g ^-x and the moduli of Xg^-x- 
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